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We examine correlation functions in the presence of competing long and short interactions to find 
multiple correlation and modulation lengths. For many systems, we show the presence of a crossover 
temperature T* above which there is no modulation and report on a new exponent, vl, characterizing 
the universal nature of this divergence of the modulation length. We study large n systems and 
find, in general, that ul = 1/2. We demonstrate that for a short range system frustrated by a 
general competing long range interaction, the crossover temperature T* veers towards the critical 
temperature of the unfrustrated short range system (i.e., that in which the frustrating long range 
interaction is removed) . We discover in systems with long range interactions the existence of at 
least one diverging correlation length in the high temperature limit. When screening is present, 
instead of diverging, this correlation length tends, at high temperatures to the screening length. 
We also show that apart from certain special crossover points, the total number of correlation and 
modulation lengths remains conserved. We derive an expression for the change in modulation length 
with temperature for a general system near the ground state with a ferromagnetic interaction and 
an opposing long range interaction. We illustrate that the correlation functions associated with the 
exact dipolar interactions differ substantially from those in which a scalar product form between 
the dipoles is assumed. 
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I. INTRODUCTION 

Short range interactions have been at the focus of much 
study for many decades. Perhaps one of the best known 
examples are the Ising ferromagnet and the antiferro- 
magnet Long range interactions are equally abun- 
dant Systems in which both long and short range 
interactions co-exist comprise very interesting systems. 
Such competing forces can lead to a wealth of interesting 
patterns - stripes, bubbles, etc [3], [1], Reahzations 
are found in numerous fields - quantum Hall systems 
a|, adatoms on metallic surfaces, amphiphilic systems 
7|, interacting elastic defects (dislocations and disclina- 
tions) in solids [8|], interactions amongst vortices in fluid 
mechanics @ and superconductors [lO|, crum pled mem- 
brane systems [llj, wave-particle interactions 12 1, inter- 
actions amongst holes in cuprate superconductors p^ . 
Hi 13 1 01 [I3|i arsenide superconductors [Tsj . man- 
ganates and nickelates Il9ll . I20l some theories of struc- 
tural gksses [U, m, m, [ij, m, conoldal systems 
[2^ . |27j and many many more. Much of the work to 
date focused on the character of the transitions in these 
systems and the subtle thermodynamics that is often ob- 
served (e.g., the equivalence between different ensembles 
in many such systems is no longer as obvious, nor al way s 
correct, as it is in the canonical short range case [2a|). 
Other very interesting aspects of different systems have 
been addressed in 29]. 

Here we investigate the general temperature depen- 
dence of the structural features that appear in such 
systems when competing interactions of short and long 
range are present. In many such systems, there are emer- 
gent modulation lengths governing the size of various do- 



mains. We find that these modulation lengths often ad- 
here to various scaling laws, sharp crossovers and diver- 
gences at various temperatures (with no associated ther- 
modynamic transition). We also find that in such sys- 
tems, correlation lengths generically evolve into modu- 
lation lengths (and vice versa) at various temperatures. 
The behavior of correlation and modulation lengths as a 
function of temperature will afford us with certain selec- 
tion rules on the possible underlying microscopic inter- 
actions. In their simplest incarnation, our central results 
as follows: 

1. In canonical systems harboring competing short (fi- 
nite) and long range interactions modulated pat- 
terns appear. Depending on the type of the long 
range interaction, the modulation length either in- 
creases or decreases from its ground state value 
as the temperature is raised. We will relate this 
change, in lattice systems, to derivatives of the 
Fourier transforms of the interactions that are 
present. 

2. There exist special crossover temperatures at which 
new correlation/modulation lengths come up or 
some cease to exist. 

(a) For many systems, there exists a sharp 
crossover temperature T* at which the mod- 
ulation length diverges. That is, we find that 
in many systems, for all temperatures larger 
than a threshold temperature, T* , there are no 
modulations in the pair correlation functions. 
If for T < T*, finite modulation lengths set in 



then these scale as 

Ld- (T* -T)-"^, (1) 

irrespective of the interaction; quite univer- 
sally, J^L = 1/2. 

(b) In some cases, the modulation length can 
change by finite increment jumps instead of 
diverging. 

(c) The total number of characteristic length- 
scales (correlation -I- modulation) remains 
conserved, except at the crossover points. 

3. In the high temperature limit of large n theories 
there generally exists one or more diverging corre- 
lation lengths, i.e., 

hm e = oo. (2) 

When screening is present, the divergence is re- 
placed by the a saturation 

hm e - X-\ (3) 

i — ^CJO 

where is the screening length. 

4. The presence of the angular dependent dipolar in- 
teraction term that frustrates an otherwise un- 
frustrated ferromagnet vis a vis a simple scalar 
product between the dipoles adds new (dominant) 
lengthscales. The angular dependence significantly 
changes the system. 

We will further investigate the ground state modula- 
tion lengths in general frustrated Ising systems and also 
point to discontinuous jumps in the modulation length's 
that may appear in the large n rendition of some systems. 

Armed with these general results, we may discern the 
viable microscopic interactions (exact or effective) which 
underlie temperature dependent patterns that are trig- 
gered by two competing interactions. Our analysis sug- 
gests the effective microscopic interactions that may drive 
non-uniform patterns such as those underlying lattice 
analogs of the systems of Fig.dT]). 

The treatment that we present in this work applies to 
lattice systems and does not account for the curvature 
of bubbles and other continuum objects. These may be 
augmented by inspecting energy functional (and their 
associated free energy extrema) of various continuum 
field morphologies under the addition of detailed domain 
wall tension forms - e.g., explicit line integrals along the 
perimeter where surface tension exists ~ and the impo- 
sition of additional constraints via Lagrange multipliers. 
We leave their analysis for future work. One of the cen- 
tral results of our work is the derivation of conditions 
relating to the increase/decrease of modulation lengths 
in lattice systems with changes in temperature. These 
conditions relate the change in the modulation length at 
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FIG. 1: Reproduced with permission from Science, Ref. [Tl[. 
Reversible "strip out" instability in magnetic and organic 
thin films. Period (Ld) reduction under the constraint of 
fixed overall composition and fixed number of domains leads 
to elongation of bubbles. Left panel (A) in magnetic garnet 
films, this is achieved by raising the temperature [labeled in 
(B) in degrees Celsius] along the symmetry axis, H — (pe- 
riod in bottom panel, ~ (see Fig. 5). Right panel (B) 
In Langmuir films composed of phospholipid dimyristolphos- 
phatidic acid (DMPA) and cholesterol (98:2 molar ratio, pH 
11), this is achieved by lowering the temperature at constant 
average molecular density [period in bottom panel, ~ 20/im] 



low temperatures to the derivatives of the Fourier trans- 
forms of the interactions present. 

In Section(|TT|, we outhne the general systems that we 
study - short range systems that are frustrated by long 
range interactions. We then proceed, in Section (jIII|l to 
introduce a crossover temperature T* and prove a gen- 
eral scaling relation associated with it. In Section pVj) . 
we derive the scaling form for the Ising ground states 
for general frustrating long range interactions. Hence- 
forth, we provide explicit expressions for the crossover 
temperatures and the correlations lengths in the large n 
limit. In Section (|V|l . we introduce the two spin corre- 
lation function for a general system in this limit. Based 
on the correlation function, we then present some general 
results for systems with competing nearest neighbor fer- 
romagnetic interaction and an arbitrary long range inter- 
action in Section (jVip . We start by deriving the equilib- 
rium stripe width for a two dimensional Ising system with 
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nearest neighbor ferromagnetic interactions and compet- 
ing long range interactions. We derive an expression for 
the change in modulation length with temperature for 
low temperatures for large n systems. We illustrate how 
the crossover temperature, T* arises in the large n limit 
and show some general properties of the system associ- 
ated with it. We show the presence of a diverging corre- 
lation length in the high temperature limit. 

We present some example systems in Section (jVIII) . 
We numerically calculate the correlation function for the 
screened Coulomb frustrated ferromagnet and the dipo- 
lar frustrated ferromagnet. We then study the screened 
Coulomb frustrated ferromagnet in more details. Next, 
we show some results for systems with higher dimen- 
sional spins. We study a system with the dipole-dipole 
interaction for three dimensional spins, without ignor- 
ing the angular dependent term and show that this term 
changes the ground state lengthscales of the system con- 
siderably. We also present a system with the Dzyaloshin- 
sky - Moriya interaction in addition to the ferromagnetic 
term and a general frustrating long range term. 

We give our concluding remarks in Section (jVIIip . 



II. THE SYSTEMS OF STUDY 

Consider a translationally invariant systems whose 
Hamiltonian is given by 



(4) 



The quantities {S{x)} portray classical scalar spins or 
fields. The sites x and y lie on a hypercubic lattice of 
size N of unit lattice constant. In what follows, v{k) and 
s{k) will denote the Fourier transforms of V^dx — y\) and 
S{x). For analytic interactions, v{k) is a function of 
(to avoid branch cuts). We will focus on systems with 
competing interactions. Throughout most of this work 
we will focus on the following Hamiltonian with short 
and long range interactions. 

H S{x)Siy)+Qj2VL{\x-y\)S{x)S{y) (5) 



where the first term represents nearest neighbor ferro- 
magnetic interaction for positive J and the second term 
represents some long range interaction which opposes the 
ferromagnetic interaction for positive Q. We will study 
properties of general systems of the form of Eq.([5]). In or- 
der to flesh out the physical meaning of our results and 
illustrate their implications and meaning, we will fur- 
ther provide explicit expressions and numerical results 
for two particular examples. The Hamiltonian of Eq.([5]) 
represents a system that we christen to be the screened 
Coulomb Frustrated Ferromagnet when 



VL{r) 



Xr 



Airr 



in three dimensions, and 



where A ^ represents the screening length and Kq is a 
modified Bessel function. 



, , f , cosxt 
Ko{x) = dt-j= 



(7) 



Throughout our work, we will discuss both the screened 
and unscreened renditions of the Coulomb frustrated fer- 
romagnet. Eq.® corresponds to a Dipolar Frustrated 
Ferromagnet when 



(^2 + ^2)3/2 



in the limit (5 — > 0, 



(8) 



on the lattices that we will consider. Later, we also con- 
sider the general direction dependent (relative to the lo- 
cation vectors) of the dipolar interaction for three dimen- 
sional spins; we will replace the scalar product form of the 
dipolar interactions in Eqs.fO E]) by the precise dipolar 
interactions between magnetic moments. 

On a hypercubic lattice, the nearest neighbor interac- 
tions in real space of Eq. ^ have the lattice Laplacian 



A(fc) = 2^(1 - coski), 



(9) 



1=1 



as their Fourier transform. In the continuum (small k) 
limit, A ^ z = fc^. The real lattice Laplacian 



{x\m = 



2d for X — y 
— 1 for \x — y\ 



1. 



(10) 



Notice that {x\A^\y) — for [x — y| > R, where R is 
the spatial range over which the interaction kernel is 
non- vanishing. Eq. (jlOp corresponds to interactions of 
Range=2 lattice constants, 

{x\A'^\y) = 2d(2d -f 2) for f = y 
—id for \x — y\ — 1 

2 for {x — y) = (±ei ± et) where £ ^ £' 
1 for a ±2ei separation. (11) 



^L{r) — Ko{Xr) in two dimensions. 



Correspondingly, in the continuum, the Lattice Lapla- 
cian and its powers attain simple forms and capture ten- 
dencies in numerous systems. Surface tension in many 
systems is captured by a g(V0)^ term where is a con- 
stant in a uniform domain. Upon Fourier transforming, 
such squared gradient terms lead to a fc^ dependence. 
The effects of curvature which are notable in many mix- 
tures and membrane systems can often be emulated by 
terms involving (V^/i) with h a variable parameterizing 
the profile; at times the interplay of such curvature terms 
with others leads, in the aftermath, to a simple short 
range k'^ term in the interaction kernel (the continuum 
version of the squared lattice Laplacian of Eq.([TT|)). An 
(6) excellent review of these issues is addressed in [llj. 



4 



III. A UNIVERSAL DOMAIN LENGTH 
EXPONENT 

We next report on a new exponent for the domain 
length of general frustrated systems. Our result below 
applies to general fields. It spans real or complex scalar 
fields, vectorial (or tensorial) fields of both the discrete 
(e.g., Potts like) and continuous variants. Although the 
result below is general, for consistency of notation in this 
work, we will refer to a system of 0{n) spins. The desig- 
nation of "0(n) spins" simply denotes real fields (spins) 
of unit length that have some arbitrary number (n) of the 
components. For n = 1, the system is an Ising model. 
A single component real field having unit norm allows 
for only two scalars at any given site x: S{x) = ±1. 
The n = 2 system corresponds to a two component spin 
system in which the spins are free to rotate in a two di- 
mensional place- Sf{x) + 51(2?") — 1 (the so-called XY 
spin system). The case of n = 3 corresponds to a system 
of three component Heisenberg type spins, and so on. In 
general, 



1 = ^ Sa{x)Sa{x). 



(12) 



In the up and coming, we will denote the momentum 
space correlator by 



G{k) = {Sik)-S{-k)).) 



(13) 



In the case of the Ising model, n — 1, the scalar prod- 
uct in Eg. p^ is to be replaced by a product between 
the complex scalars S{k) and S{—k). The characteristic 
lengthscales of the system are given by the solutions of 
G~^{T, k) — where G~^{T, k) denotes the inverse pair 
correlator at wavenumber k and temperature T. We will 
now consider the general case of Eq.([3]) in which the sys- 
tem exhibits modulations at a commensurate wavenum- 
ber q for all temperatures T > T* , and starts to exhibit 
incommensurate modulations at temperatures T = T*^ . 
Throughout this work, we will provide several specific 
examples. Whenever the above conditions are satisfied 
(i.e., incommensurate modulations set in below a tem- 
perature T*) then the modulation must universally follow 
the power law scaling as a function of the temperature 
(T* — T) with an exponent of vl = 1/2. We first show 
how this result occurs generally. 

As the system exhibit modulations of wavenumber q 
above T*, the real part of the poles of the correlator 
must be q. That is. 



G-^{T,q + iK) ^ 0. 



(14) 



If plot the curve for Eq.lHl]) in the T - k plane, T* will 
denote the point where the curve just touches the line 
T = T* from above. The curve is parallel to the K-axis 
at that point. Thus, 



QJ^ \(T*.q+tK) 



0. 



(15) 



We investigate the case in which T* ^ and for which 
the inverse correlator is analytic at T* and an expansion 
is possible. For poles at a slightly lower temperature, we 
have to lowest order. 



G -\T,k)^0 



{T* - T) 



dG-^ 
dT 



\{T',q+iK) 



-\{T',q+iK) 



= 0. 



2! 9fc2 

Thus, if Sk is not imaginary then, to leading order. 



(16) 



(X {T* -TY^, 



with 



(17) 



Eq. ()17|) constitutes a new counterpart for the well known 
divergence of the correlation length at Tc. There is now 
a new divergence of the modulation length at the "disor- 
der" line temperature T* . [s^ The value of the critical 
exponent is similar to that appearing for the correlation 
length exponent for mean-field or large n theories. It 
should be stressed that our result of Eq. p7|) is far more 
general. 



IV. GROUND STATE STRIPE WIDTH FOR 
ISING SYSTEMS 

We next briefly address the opposite extreme - low 
temperatures, and discuss the ground states of general 
long range interactions in Eq.Q. Below, we discuss the 
Ising ground states. We will later on consider the spher- 
ical model that will enable us to compute the correlation 
functions at arbitrary temperatures. We consider a sys- 
tem with Ising spins in d dimension and assume that the 
system forms a "striped" pattern (periodic pattern along 
one of the dimensions - stripes in two dimensions, planes 
that lie in three dimensions and so on) of spin-up and 
spin-down states of period /. We then calculate the free 
energy as a function of I and then minimize it with re- 
spect to I to get the equilibrium stripe width. For the 
frustrated ferromagnet, if VL{k) — l/k^, 



I 



(27r)P+V 



(18) 



where ("(5) is the Riemann zeta function, 

00 

C(.s) = ^n-. (19) 

n=l 

For Coulomb frustration, p = 2, 



I = 



'647r2 J 



(20) 
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in accord with the results of Refs. [23, l3l|. For long 
range dipolar interactions {p ~ 3), we find that 



(21) 



In the notation to be employed throughout this work, I 
plays the role of the modulation length of the system at 
zero temperature, Lu{T — 0). 



V. CORRELATION FUNCTIONS IN THE 
LARGE N LIMIT - GENERAL 
CONSIDERATIONS 

The results reported henceforth were computed within 
the spherical or large n limit . It was found by Stanley 
long ago that the large n limit of the n component nor- 
malized spin systems (so called 0{n) spins) introduced 
in Section piI[) is identical to the spherical model first 
introduced by Berlin and Kac. |32] We now introduce 
this model in its generality. The spins in Eq. Q satisfy a 
single global ("spherical") constraint, 



(22) 



enforced by a Lagrange multiplier /i. This leads to 
the functional H' = H + jiN which renders the model 
quadratic (as both Eqs.(l4l [22l) are quadratic) and thus 
exactly solvable, see, e.g., [14| . The continuum analogs 
of Eqs.(|l22]) read 

H=^Jd''x d'^y V{\x - y\)S{x)S{y), 

d'^x {S^{x)) ^ const. (23) 



The quadratic theory may be solved exactly. From the 
equipartition theorem, for T > Tc, the Fourier space cor- 
relator 



G(fc) ^ {\S'{k)\) 



v{k) + /i 

The real space two point correlator is given by 



G{x) = (S'(O)S'(f)) = keT 



d<^k 



BZ 



{2itY v{k)+n' 



(24) 



(25) 



with d the spatial dimension and BZ denoting the inte- 
gration over the first Brillouin zone. For a hypercubic 
lattice in d dimensions with a lattice constant that is set 
to one, — TT < ki < TT for i ~ l,2,...,d. Henceforth, to 
avoid cumbersome notation, we will generally drop the 
designation of BZ; this is to be understood on all mo- 
mentum space integrals pertaining to the lattice systems 
that we examine. To complete the characterization of the 
correlation functions at different temperatures, we note 



that the Lagrange multiplier /i(T) is given by the implicit 
equation 1 ~ G{x — 0). Thus, 



l = kBT 



d'^k 



1 



(27r)''i;(fc) + ^i' 



(26) 



This implies that the temperature T is a monotonic in- 
creasing function of jj.. If fj, changes by a small amount 
A/x, then T will change by an amount AT, such that 



Eq. 



AT cx A^. (27) 
implies that in the high temperature limit, 



kBT 



Bi (2^)'*r(f 
T cx /X, 



1) 



(28) 
(29) 



where A represents the upper limit of the k integration, 
representing the ultra-violet cut-off. If we perform the 
momentum integration in a hypercube of side 27r, we have 
in the high temperature limit. 



^ = kBT. 



(30) 



We investigate the general character of the correlation 
functions given by Eq. (|25p for rotationally invariant sys- 
tems. If the minimum (minima) of v{k) occur(s) at mo- 
menta q far from the Brillouin zone boundaries of the 
cubic lattice then we may set the range of integration 
in Eq. (P5)) to be unrestricted. The correlation function 
is then dominated by the location of the poles (and/or 
branch cuts) of l/[v{k) + /i]. Thus, we look for solutions 
to the following equation. 



v{k) + /i = 0. 



(31) 



The system is perfectly ordered in its ground state. 
From a temperature at which the system is not perfectly 
ordered, as we lower the temperature, the correlation 
length diverges aX T — Tc- For T = Tc, fj, takes the 
value. 



^J.nlin = - min [v{k)]. 

keBZ 



(32) 



As the temperature is increased, the disorder creeps in 
and in most systems, at a temperature T*, the modula- 
tion length diverges. 

The characteristic length scales of the system are gov- 
erned by the poles of [v{k) + ^]~^. 

JA{k) + QvL{k)+ n^O, (33) 

which in the continuum limit takes the form 

Jk^ + QvL{k) + ^ 0. (34) 

Employing the above considerations, we will derive, in 
the next section, some general results for systems of the 
form Eq.®. 
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Our work will focus on classical systems. The exten- 
sion to the quantum arena [T^ is straightforward. In, 
e.g., large n bosonic renditions of our system, replicat- 
ing the usual large n analysis, [1^ we find that the pair 
correlator 



v{k)+ii 
ksT 



with the bosonic distribution function 



nBi-j-^) = [exp(— — ) - 1] 



keT' 



(35) 



(36) 



The correlator of Eg. ([35)1 is of a similar nature as that 
of the classical correlator of Ea. ([M|) with branch cuts 
generally appearing in the quantum case. Our analysis 
below relies on the evolution of the poles of v{k) -\- jj. 
as a function of temperature in classical systems. In the 
quantum arena, these poles are replaced by branch points 
and the analysis remains qualitatively the same. 



VI. RESULTS 

In this section, we present some general results for sys- 
tems of the form Eq.([5]) in their large n limit. First, 
we find the dependence of the modulation length on 
temperature, near Tc- Next, we will illustrate an anal- 
ogy between the behavior of the correlation length near 
the critical temperature and that of the modulation 
length near T* . We will then discuss some aspects of the 
crossover points. We end this section with some results 
in the high temperature limit. 



A. The low temperature limit: a criterion for 
determining an increase or decrease of the 
modulation length at low temperatures 

In this section, we derive universal conditions for in- 
creasing or decreasing modulation lengths in systems 
with competing ferromagnetic and long range interac- 
tions. Eqs. (|T7[B5)l show general conditions for the change 
in modulation length, Ljj with temperature for a general 
system of the form Eq.Q. The value, fep of k which 
satisfies Ea. (|32)) . 



w(fco) ~ min vik) 

keBZ 

determines the modulation length at T = Tc- 

Jkl+ QvL{ka) + ^imin = 0, 

2Jko + Qv'i^{ko) =0. 



(37) 



(38) 
(39) 



As the temperature is raised, the new pole near fco will 
have an imaginary part corresponding to the finite cor- 
relation length. The real part will also change in general 



and this would induce a change in the modulation length. 
Let fi = fJ-min + Sfi with Sfj, > 0. Then wc have. 



k ~ ko + Sk, 

oo 

Sk = 6k j, 



(40) 



where Skj cx ^/i^^, Xj+i > Xj. Our goal is to find the 
leading order real contribution to 5k which would give us 
the change in modulation length with increasing Sfi and 
hence with increasing temperature. 



JSk' + Qf2v^^\ko)^-^+S^l = 0. 



(41) 



Suppose ^^^"■'(fco) = for 2 < n < p and v'f ^ko) =/= 0. 
We have, 

(J+^^^sJ^){Skl + 2Sk,Sk, + ...) 



2! 

'^''''^^''"ksk{+p6kr'6k. 



p 



ip+iy- 



iSkP+^ + {p+ l)SklSk2 + ...) + ...] +6^ = 0. (42) 
To leading order. 



{J 



2! 

6ki 



6^ 



(43) 



J 



From this, we see that 5ki is imaginary. This consti- 
tutes another verification of the well established result 
about the universality of the divergence of the correlation 
length, ^ at Tc with the mean-field type critical exponent 
V — 1/2 in the large n limit. 



C(x(^-^,)-^ 



(44) 



The next, higher order, relations are obtained using the 
method of dominant balance. 

Therefore, 5k2 is real if p is odd and imaginary if p is 
even. If, 



Ld{T)=.Ld{Tc) + 6Ld, 
then, for p = 2n + 1 , 

2n (-l)"Q4^"+^>(fco)V 



2(2n + l)!(J+^^;^-y^)»+i 



(46) 



(47) 
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Thus to get the leading order real contribution to 5k for 
even p[> 2], we have to go to higher order. 



{-ly+p/^Qv^p+'\ko){S^l)p/' 



For p = 2n, 



SLj 



2tt 



2(p+l)!(J + 



(-1)"Q4'"^ 



2! 



)p/ 



2+1 



\ko){5^xy 



^0 2(2n+l)!(J+^ 



(48) 



(49) 



n+l 



If, for p = 2n, w^^"^^'(fco) = 0, then we will need to 
continue this process until we get a real contribution to 
5k. In appendix |Bl we provide explicit forms for 5L]j 
for different values of p. The results derived above al- 
low us to relate an increase /decrease in the modulation 
length at low temperatures to the sign of the first non- 
vanishing derivative ( of an order larger than two ) of the 
Fourier transform of the longer range interactions that 
are present. It is important to emphasize that our re- 
sults apply to any interaction, vl that augments a nearest 
neighbor type interaction. These may include screened 
or unscreened Coulomb and other long range interactions 
but may also include interactions that are strictly of finite 
range (e.g., next-nearest neighbor interactions on the lat- 
tice for which we have vl = —tA^ (with a constant i > 0, 
see Eq.dlll))). 

The results from this section about modulation 
lengths, can give us similar behavior of the correlation 
lengths at temperatures slightly below T*. 



B. A correspondence between the temperature T* 
at which the modulation length diverges and the 
critical temperature Tc 

The critical temperature Tc corresponds to the max- 
imum value of fi for which Eq.(|3ip still attains a real 
solution. Thus, 

v{ka) + ^Irnin = 0, 

f'(fco) =0, 

v'ika) > 0. (50) 

For systems in which the modulation length diverges at 
T* , T* corresponds to the minimum value of /i for which 
Eq. ()31|) has a purely imaginary solution. Thus, if v{iK) = 
v{k), 



v"inQ) < 



v{ko) + /i* = 
i)'(Ko) = 
=> v"{iKo) > 



(51) 



Thus, we expect similar qualitative results for the correla- 
tion lengths at temperatures slightly above as for mod- 
ulation lengths slightly below T* and vice-versa. [The 



relations for the derivatives of v{ko) in Eq. dSTI) are guar- 
anteed to hold only if T* > T^.] 



C. Crossover temperatures: Emergent modulations 

For systems with competing multiple range interac- 
tions, there may exist special temperatures at which the 
poles of the correlation function change character, thus 
changing modulation lengths to correlation lengths and 
vice-versa. In particular, for most systems we have a 
crossover temperature T* above which the system does 
not have any modulation. Apart from this kind of phe- 
nomenon, there might also be finite discontinuous jumps 
in the modulation length. This is illustrated with an ex- 
ample in Section (jVIEp . 

We start by defining the crossover temperature T* for 
a ferromagnetic system frustrated by a general long range 
interaction. Let k = iK, k £ R above T* and k = kq at 
T*. Let v{k) = f{z), z = fc^. Above T* , ^ = /i„i„ -f A/i 
[A^ > 0]. Using Eq.JSll), 



u = — f(—K ) = — minw(A;) - 
-v{k)] 



max 

km 



f A/1. 



(52) 



T* corresponds to the minimum value of A/i for which 
we have at least one such solution. See Fig. ([2]). Thus, 



50 



25 



^ 
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mm 











-4 -3 -2 -1 



FIG. 2: —f{z) — —v{k) = fi plotted against z = for 
purely real and purely imaginary fc's (T — > and T — > oo 
respectively). The negative z regime corresponds to temper- 
atures (Lagrange multiplier, /i's) for which purely imaginary 
solutions exist. The maximum of the curve in the positive z 
regime corresponds to the modulations a.t T = Tc [^i = /imm], 
which is the maximum temperature at which pure modula- 
tions exist. 



mm 
— max 

— v{iK) > IJr, 



-v{iK)] 
[viin)]. 



(53) 
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Sometimes, the crossover point is slightly more difficult 
to visualize. See Fig. ([3]). In this case, the minimum up- 
per branch of —f{z) for z < [equivalently the upper 
branches of —v{iK)] gives us the value of /i*. The branch 
chosen has to continue to ii — +oo so that at least one 
term without modulation is always available as we in- 
crease the temperature, as required by the definition of 
T*. The other branch provides such solutions only up to 
a certain temperature. Also, the part of it which is below 
^J■m,^n IS irrelevant. 



— f{z) to exist in the z < regime which is higher than 
the maximum in the z > regime, we need f{z) to be 
discontinuous at some point. Putting all of the pieces 
together, we find that there are no possibilities: (i) no 
crossover, i.e., T* = oo or (ii) kq = and /i* — fimin 
with T* = Tc. 



2. T* Tc as the strength of the long range interaction is 
turned off 
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FIG. 3: TOP: Solid line:~-!;(fe) plotted against k. 
Dashed line:— plotted against k. 
BOTTOM: -f{z) plotted against z. 



If f{z) is an odd function of z (e.g. the Coulomb frus- 
trated ferromagnet), /i* = — /imm and the correlation 
length at T* is the same as the modulation length at Tc- 

Also, for the system in Eq.([5]), if limK_,o VLiin) — -f cxd, 
we have, /i* = ^imin, and T* — T^. 



The results from this section and the next hold for a 
general system, not just the frustrated ferromagnet. 

The crossover temperature T* tends to Tc for Q = 
as Q ^ 0. For a general system, let G(T, k) denote the 
Fourier space correlation function at temperature T. By 
definition, at T = Tc the correlation length is infinity. 
Thus, Tc is the solution to 



G-\T,k) = Q, 



(54) 



such that k G BZ (or for continuum renditions, fc G K). 

T* is the temperature at which the modulation length 
diverges for the frustrated ferromagnet, or becomes the 
same as the modulation length of the unfrustrated system 
at Tc for a general system. Thus, T* is the solution to 



G-HT,q + iK) = 0, 



(55) 



with K G M [(7 = for the case of the frustrated ferromag- 
net, q ~ t: for the frustrated antiferromagnet] . At Tc, for 
Q — 0, we have. 



G-\Tc,q)=0. 



(56) 



This however also satisfies Eg. (1551) . with k = 0. There- 
fore, 



lim T* =Tc. 



(57) 



We demonstrate this in the large n limit. See Fig.([l]). 
For Q = 0, we have ^J-miri = and fco = 0. Let VL{k) 
diverge as k~^P near fc = 0. For small Q, from Eg. ([5^ . 
we have. 



fJ-mzn = -^Jp+iQp+i. 



(58) 



1. T* = Tc if all the competing interactions are of finite 
range and crossover exists 

For systems where all the competing interactions are 
of finite range, T* = Tc- We prove this as follows. Since 
finite range interactions contribute to v{k) as powers of 
A(fc) — !■ fc^, for a system with only finite range inter- 
actions, f{z) is analytic for all z. For a minimum of 



If p is odd, 



(59) 



As Q — > 0, = fcg = and /i* = l-tmin = Oi that is, 

lim T* ^TciQ^Q). (60) 



9 



3. Proof of the conservation of the total number of 
characteristic lengthscales 

In this section we consider the general situation in 
which the Fourier transform of the interaction kernel, 
v{k), is a rational function oi z — A(fc), [z — > fc^ in 
the continuum limit]. That is, we consider situations for 
which 

^(fc) = /W = ^, (61) 

with P and Q polynomials (of degrees Mi and M2 re- 
spectively). We will now demonstrate that the combined 
sum of the total number of correlation and the number 
of modulation lengths remains unchanged as the temper- 
ature is varied. 

Before providing the proof of our assertion, we first re- 
iterate that the form of Eq. (pTj) is rather general. For a 
system with finite range interactions (t^da'— yl > i?) = 
with finite R) that is even under parity {V{x—y) = V{y— 
x)), the Fourier transform of V{x — y) can be written as 
a finite order polynomial in (1 — cos/c/) with the spatial 
direction index 1 < I < d where d is the dimensionality. 
In the small \k\ (continuum limit), [1 — cosfc;] kf/2. 
The particular case of a system with only finite range 
interactions that exist up to a specified range R on the 
lattice (the range being equal to a graph distance mea- 
suring the number of lattice steps beyond which the in- 
teractions vanish) of the form of Eq. (pTjl corresponds to 
v{k) = P{z) with the order of the polynomial Mi being 
equal to the interaction range, R = Mi. Our result be- 
low includes such systems as well as general systems with 
long range interactions. For long range interactions such 
as, e.g., the screened Coulomb frustrated ferromagnet, 
f{z) = 1/(2 + A^). The considerations given below apply 
to the correlations along any of the spatial directions I 



-2-1 1 2 



z 

FIG. 4: Illustration of the limit T* ^ Tc a.s Q ^ with 
VL{k) — l/k'^. The plot shows —f{z) — —v{k) vs z = , for 
v{k) = Jk^+Q/k^ with J = 1 and Q = {1-Blue, 0.1-Green, 
0.01— Yellow, 0.001— Red}. represents the value of p* and 
dot represents Umin- 



(and as a particular case, radially symmetric interactions 
for which the correlations along all directions attain the 
same form). 

Returning to the form of Eq. (pT|) . the Fourier space 
correlator of Eq. ((24|) is given by 

G(^)-fcsT||^; F{z)^P{z)+^,Q{z). (62) 

On Fourier transforming Eq. ([5^ to real space to ob- 
tain the correlation and modulation lengths, we see that 
the zeros of F{z) determine these lengths. Expressed in 
terms of its zeros, F can be written as 

M 

F{z)=A\{{z-z,), (63) 

where M = max[Mi,M2]. Perusing Ea. ([S^ . we see 
that F is a polynomial in z with real coefficients. As 
F*[z) = F{z*) it follows that all roots of F are ei- 
ther (a) real or (b) come in complex conjugate pairs 
{zj = z* ^ Zi). We now focus on the two cases sepa- 
rately. 

(a) Real roots: If a particular root = a^, a G M then 
on Fourier transforming Eq. ([62)1 by the use of the residue 
theorem, we obtain a term with a modulation length, 
Lu = l-nja. Conversely, if Zj = —a?, we get a term with 
a correlation length, <^ — 1/a. 

(b) Next we turn to the case of complex conjugate 
pairs of roots. If the pair of roots Zj,z* is not real, 
that is, Zj — |zj|e'^, then on Fourier transforming, 
we obtain a term containing both a correlation length, 
^ = (\/|zj|| sin ||)~^, and modulation length, Ld = 
2^(v^|cosf|)-i. 

Putting all of the pieces together see that as (a) each real 
root of F{z) contributes to either a correlation length or 
a modulation length and (b) complex conjugate pairs of 
non-real roots contribute to one correlation length and 
one modulation length, t/ie total number of correlation 
and modulation lengths remains unchanged as the tem- 
perature (fj.) is varied. The total number of correlation 
-|- modulation lengths is given by the number of roots 
of F{z) (that is, M). Thus, the system generally dis- 
plays a net of M correlation and modulation lengths. 
This concludes our proof. At very special temperatures, 
the Lagrange multiplier ^(T) may be such that several 
poles degenerate into one - thus lowering the number of 
correlation/modulation lengths at those special temper- 
atures. Also, in case M = M2, the total number of roots 
drops from M2 to Mi at fi = 0. What underlies mul- 
tiple length scales is the existence of terms of different 
ranges (different powers of z in the illustration above) - 
not frustration. 

The same result can be proven using the transfer ma- 
trix method, for a one-dimensional system with Ising 
spins. This is outlined in appendix |Al A trivial exten- 
sion enables similar results for other discrete spin systems 
(e.g.. Potts spins). 
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D. The high temperature limit 

In this section, we show that in the high temperature 
hmit, in systems with competing short and long range in- 
teractions, there exists the usual correlation length which 
goes to zero. Apart from this, we show that there also 
exists at least another correlation length which diverges 
or tends to the screening length if screening is present. 
The correlation function however satisfies, 



lim G{x 7^ 0) = 0. 



(64) 



Thermodynamically, the pre-factor multiplying a term 
with diverging correlation length should go to zero in the 
high temperature limit and it indeed does. 

The results in this section are valid for an arbitrary 
long range frustrating interaction VL{r), in the large n 
limit. 

The correlation and modulation lengths at the high 
temperature is governed by the roots of Ea. (p4)l . for large 
fi. The long range interaction usually satisfies, 

lim tiL(/c) = oo, (65) 
lim WL(fc) = 0. (66) 



We look for self-consistent solutions to Ea. (p4| 



1. Small correlation length ~ The usual behavior 

Small correlation lengths correspond to a large |fc| 
of the poles of the pair correlator in Fourier space (or 
more precisely, to a large imaginary part of the poles 
in k space). We now look for self-consistent solution to 
the equation G^^(fc) = that have large |fc| and for 
which |fcp ^ l''^L(fc)l £ind because the temperature is 
high and the system has a high value of fi also satisfy 
A* 3> |wL(fc)|- With all of the above in tow, the pole equa- 
tion G~^{k) = reduces to the equation obtained in the 
"canonical" case of short range interactions. 




(67) 



Thus, at least one correlation length is, to leading order. 



2. A new divergence or saturation (in presence of 
screening) of the correlation length 

We search for a self-consistent solution with a small 
modulus |fc| <C 1 to the equation G^^(fc) = in the 
limit of high temperatures (high /i). A small value of |fc| 
mandates a small value of the imaginary part of k and 
thus to a large correlation length. Setting |fL(fc)| ^ 
(and /i ^ |/cp) wc find, to leading order, that G^^(fc) = 
reads 



QvL{k) = -fi, 



(70) 



where, 
relation length, 



is the inverse function of v^. That is, 
= q. From Eq. ([65|) . we find that the cor- 



(71) 



We next briefly comment on the momentum space kernel 



VLik) 



(fc2 + A-2)p/2 



(72) 



that corresponds to a screened variant of V{r) ^ ^ with 
X < d for p > with a screening length A^^ which we 
take to be small. In the high temperature limit, we find 
a self-consistent small solution to the pole equation 
G~^{k) = 0. This solution illustrates that 



lim 1^ = A 



(73) 



That is, at least one of the correlation lengths tends to 
the screening length A~^ in the limit of high tempera- 
tures. In the absence of screening, i.e., A~^ = 0, at least 
one correlation length diverges in the high temperature 
limit! This might seem paradoxical. However, a careful 
evaluation shows that in this limit the prefactor associ- 
ated with this correlation length tends to zero as T ^ oo 
and all correlations decay monotonically with increasing 
temperature as they must. The reader can see how this 
occurs in some of the example systems that we will ex- 
amine in detail later on (e.g., Eqs. ((85l [861) 1. 



E. First order transitions in the modulation length 



(68) 



Ea. (j29p states that in the high temperature limit, in any 
dimension, 



lim £ cx — 

T-oo Jt 



(69) 



Such a behavior is anticipated; at high temperatures, the 
exhibits correlations only up to small distances. 



In this section, we show that there might be systems in 
which the modulation length makes finite discontinuous 
jumps. In these situations, the modulation length does 
not diverge at a temperature T* (or set of such tem- 
peratures). The ground state modulation lengths (the 
reciprocals of Fourier modes {iji} minimizing the inter- 
action kernel) need not be continuous as a function of 
the parameters that define the interactions. As we will 
simply illustrate below, in a manner that is mathemati- 
cally similar to that appearing in the Ginzburg-Landau 
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constructs, a "first order transition" in the value of the 
ground state modulation lengths can arise. Such a pos- 
sibility is quite obvious and need not be expanded upon 
in depth. As an illustrative example, let us consider the 
Rangers interaction kernel 



v{k) = a[A 



(74) 



with [0 < e <C 1] and c > 0. If the parameters are 
such that a > and 6 < 0, then v{k) displays three 
minima, i.e. [A + e] =0 and [A + e] = ±to^ , where 

= ■^[—b + \/\P- — 4ac]. the locus of points in the 
ah plane where the three minima are equal to one an- 
other is defined by w(fc) = 0. This leads to the rela- 
tion TO^ — ~X' Putting all of the pieces together, we 
see that h — —A^Jca/2> constitutes a line of "first or- 
der transitions". On traversing this line of "first order 
transitions" , the minimizing [A -I- e] (and thus the mini- 
mizing wavenumbers) changes discontinuously by Am = 



VII. EXAMPLE SYSTEMS 

In this section, we will investigate in detail several 
frustrated systems. We will start our analysis by ex- 
amining the screened Coulomb Frustrated Ferromagnet. 
A screened Coulomb interaction of screening length 
has the continuum Fourier transformed interaction ker- 
nel v{k) = [fc^ + A^]~^. The lowest order non- vanishing 
derivative of v^ik) of order higher than two is that of 
p — i. Invoking Eq. (j47p . we find a modulation length 
that increases with decreasing temperature as T ^ T+ 
(see also appendix B, Eq. (jBl[) in particular) . 

The dipolar interaction can be thought of as the (5 — > 
limit of the interaction, 



1 



[(x-y)2+,52]3/2- 



(75) 



This form has a simple Fourier transform. In two spatial 
dimensions, 



-kS 



(76) 



In three dimensions. 



Vd{k) = AwKoikS), 



(77) 

In this case, we similarly find that the first non- 
vanishing derivative of vl is order of order p = 3 in the 
notation of Eq. (j47| . This, as well as the detailed form 
of Eq. (|B1|) suggest an increasing modulation length with 
decreasing temperature as T ^ T+. 



A. Numerical evaluation of the Correlation 
function 

In Figs. (|5l6p . we display a numerical evaluation of 
the correlation function for the Coulomb frustrated fer- 



romagnet and the dipolar frustrated ferromagnet (see 
Eqs.([5l[ni[H])) in a two dimensional region of size 100x100. 
In both systems, we see that for low temperatures, the 
modulation length increases with temperature, thus sat- 
isfying our criterion of Eq. ((47|) [see also appendix B, Eq. 
(|B1|) in particular]. 



B. Coexisting short range and screened Coulomb 
interactions 

In this section, we study the screened Coulomb frus- 
trated ferromagnet in more details. The Fourier trans- 
form of the interaction kernel of Eq. ^ is 

Q 



j{k) = jk' + 



k^ + X^ 



(78) 



In Appendix (C), we provide explicit expressions for the 
dependence of fi on the temperature T. This dependence 
delineates the different temperature regimes. For T > T* 
wherein the temperature T* is set by fJ,{T*) = JX^ + 
2^/JQ), the pair correlator in c? = 3 dimensions is given 
by 

G(f) 



1 



47rJ|f| 



X e 



l(A2 



-I3\x\ 



Here, 



2 a2 



A2 + fi/J T ^iX^-fi/jy~4Q/J 



(79) 



(80) 



By contrast, for temperatures T < T*, we obtain an 
analytic continuation of Ea. (|79p to complex a and f3, 



G{x) 



ksT 



-ai\x\ 



8aia2TTJ\x\ 

X [(A^ — al + o?r^ sina2|a;| -f 2a\oi2 cosa2|a;|], . (81) 

In Eq. (j8ip . a — oi\ ia2 — [3* . In a similar fashion, in 
d = 2 spatial dimensions, for T > T*, 

keT 



G{x) = 



1 



27r /32 



~a^)Ka{a\x\) 
A2)Ko(/3|x|)]. 



(82) 



As in the three dimensional case, the high temperature 
correlator of Eg. ([82]) may be analytically continued to 
lower temperatures, T <T*, for which a and (3 become 
complex. 



High temperature limit 

In the high temperature limit, in two spatial dimen- 
sions. 



2ttJ' 



AttJ 



ksTA'i 



(83) 
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FIG. 5: The correlator G{x,y) for a two dimensional Screened Coulomb Ferromagnet of size 100 by 100. J — 1, Q = 0.0004, 
Screening lengtli= 100\/2. A: n = Umir. = -0.0874, B: fi = fi^in + 0.001, C: ^ = fj.min + 0.003, D: G(x,y) for j/ = for 
A(blue)[LD = 20], B(green)[Li3 = 24] and C(red)[LD = 26]. 



In three spatial dimensions, 



G{x) 



-\\x\ 



47rJ|f| kBTA^\x\ 
In the unscreened case, in two spatial dimensions, 



(84) 



G{x) 



2ttJ 



In three spatial dimensions, 
r ksT 

9tt^Q 
kBTA% 




(85) 



JiiiQ 1^1 



(86) 



From the above expressions, it is clear that the coeffi- 
cients of the terms corresponding to the diverging corre- 
lation length goes to zero in the high temperature limit. 

We note that two correlation lengths are manifest for 
all (/^ — JX^Y > ^JQ- This includes all unfrustrated 



screened attractive Coulomb ferromagnets (those with 
Q < 0)). The evolution of the correlation functions may 
be traced by examining the dynamics of the poles in the 
complex k plane as a function of temperature. At high 
temperatures, correlations are borne by poles that lie on 
the imaginary k axis. 



Thermal evolution of modulation length at low temperatures 



At T = T*, the poles merge in pairs at k 



±i\J X"^ + \JQIJ- At lower temperatures, T < T* , the 
poles move off the imaginary axis (leading in turn to os- 
cillations in the correlation functions). The norm of the 
poles, jaj = {Q / J + X"^ ii{T) / jy/^ tends to a constant in 
the limit of vanishing screening (A^^ = 0) wherein the 
after merging at T = T*, the poles slide along a circle 
[Fig. ([7])]. In the low temperature limit of the unscreened 
Coulomb ferromagnet, the poles hit the real axis at fi- 
nite fc, reflecting oscillatory modulations in the ground 
state. In the presence of screening, the pole trajectories 
are slightly skewed [Fig.®] yet for Q/ J > A^, a tends to 



the ground state modulation wavenumber y \/Q/J — X^ . 
If the screening is sufficiently large, i.e., if the screening 
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length is shorter than the natural period favored by a 
balance between the unscreened Coulomb interaction and 
the nearest neighbor attraction (A > {Q/ J)^/^), then the 
correlation functions never exhibit oscillations. In such 
instances, the poles continuously stay on the imaginary 
axis and, at low temperatures, one pair of poles veers to- 
wards fc = reflecting the uniform ground state of the 
heavily screened system. 

To summarize, at high temperatures the pair corre- 
lator G{x) is a sum of two decaying exponentials (one 
of which has a correlation length which diverges in the 
high temperature limit). For T < T* in under-screened 
systems, one of the correlation lengths turns into a modu- 
lation length characterizing low temperature oscillations. 
At the cross-over temperature T*, the modulation length 
is infinite. As the temperature is progressively lowered, 
the modulation length decreases in size - until it reaches 
its ground state value. The temperature T*{Q/J,X) is 
a "disorder line" like temperature [Fig([9])]. An ana- 
lytical thermodynamic crossover does occur at T = T*. 
A large n calculation illustrates that the internal energy 
per particle 

^ = l(fcsT-M), (87) 
To detect a crossover in U and that in other thermody- 



namic functions, the forms of both above and below T* 
may be derived from the spherical model normalization 
condition to find that the real valued functional form of 
^(T) changes [See appendix [C] . 

The system starts to exhibit order at the critical tem- 
perature T = Tc given by 

1 _ f 1 
ksTc J {2^Y v{k) - v{q)' 

For Q/J > A*, the modulus of the minimizing (ground 
state) wavcnuniber (|q|) is given by 

1=^ = ^JVQJJ-^^. (89) 

with L|) the ground state modulation length. Associated 
with this wavenumber is the kernel v{q) — — JA^ 

to be inserted in Eq. (|88|) for an evaluation of the critical 
temperature Tc- Similarly, the ground state wavenumber 
q = whenever Q/J<X^. li Q/J > and modu- 
lations transpire for temperatures T < T*, the critical 
temperature at which the chemical potential of Eq. ((25)) . 
/i(Tc) — JA^ — 2y/JQ, is lower than the crossover tem- 
perature T* (given by /i(r*) = JA^ + at which 
modulations first start to appear. The Screened Coulomb 
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T<T T <T<T* T=T* T>T* 



FIG. 7: Location of the poles with increasing temperature (left to right) in the complex k-plane for the Coulomb frustrated 
ferromagnet. For temperatures below Tc, all the poles are real. Above Tc, the poles split in opposite directions of the real 
axis to give rise to two new complex poles. For Tc < T < T* , we have complex poles. At T* , pairs of such poles meet on the 
imaginary axis. Above T* , the poles split along the imaginary axis. Thus, above T* , the poles are purely imaginary. 



k k k k. 




1=0 o<;i<(Q/j)^^'' A,= (Q/j)^^'' xxQ/j)^^'^ 

FIG. 8: Trajectory of the poles in the complex fc-plane for Tc < T < T* for the screened Coulomb ferromagnet. The screening 
length, decreases from left to right. In the first figure A = and A > (Q/J)^^'^ in the last figure. 



Ferromagnet has T^Q/J = A'') > in d > 4 dimensions 
and in any dimension T^^Q/J > A^) = 0. For small fi- 
nite n, a first order Brazovskii transition may replace the 
continuous transition occurring at Tc within the large n 
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FIG. 9: Temperature at which the modulation length diverges 
for a 100 X 100 Coulomb frustrated ferromagnet plotted versus 
the relative strength of the Coulomb interaction with respect 
to the ferromagnetic interaction. [Blue:A = Ao = l/(100\/2); 
Red:A = 0.999Ao; Black:A = l.OOlAo] 



limit 34]. Depending on parameter values such an equi- 
librium transition may or may not transpire before a pos- 
sible glass transition may occur (2^ . 

Domain length scaling in the Coulomb Frustrated 
Ferromagnet 

The characteristic length scales are governed by the 
position of the poles of [v{k) + ■ See Fig.® for an 
illustration of the pole locations at low temperatures. For 
the frustrated Coulomb ferromagnet of Eq. (|75|) in the 
absence of screening (A~^ = 0), 

v{k) + ^^=^{k^ + l^k' + ^). (90) 

Ea. (j90p enable us to determine, in our large n analysis, 
the cross-over temperature T* at which jj,* — fi(T*) — 
2\J JQ. At T = T* , the poles lie on the imaginary axis in 
fc-plane. As the temperature is lowered below T*, the two 
poles bifurcate. This bifurcation gives rise to finite size 
spatial modulations. At temperatures T < T* , the four 
poles slide along a circle of fixed radius of size (Q/ JY^'^ 
(see Figl?]). At zero temperature, these four poles merge 
in pairs to form two poles that lie on the real axis. The 
inverse modulation length is set by the absolute values of 
the real parts of the poles. We will set fi = {2y/JQ — Sfi). 
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In the following, we will obtain the dependence of the 
real part of the poles on (5/i. The poles of l/{v{k) + ji) 
are determined by 



'^pole 




±2ie 



(91) 



At /X = /i*, the angle 9 — tt/2. This point corresponds 
to the transition between (i) the high temperature region 
[T > T*) wherein the system does not exhibit any mod- 
ulations and (ii) the low temperature region (T < T*). 
[See Fig.®.] Eq.dlll) implies that cos 261 = [1 - ^] or 



kpole 



2J1/2' 

in accord with the general result of Eq. ([TT]) . 



(92) 



C. Full direction and location dependent 
dipole-dipole interactions 

In this subsection and the next, we consider systems 
where the spins are three dimensional and the interac- 
tions have the appropriate directional dependence. In 
this subsection, we will consider the effect of including 
the full dipolar interactions vis a vis the more commonly 
used scalar product form between two dipoles that is per- 
tinent to two dimensional realizations. The dipolar in- 
teraction is given by 



Hdtp - ^ 

S^y L 



S{x)-S{y) 3[S{x) ■ r\[S{y) ■ r\ 



(93) 



The two point correlator for a ferromagnetic system frus- 
trated by this interaction is given, in the large n approx- 
imation, by 



G{x) ^ keT 



(2^ 



J/\{k)^Qvd{k)+ 11 
1 

JA(fc) - 2Qvd{k) + n 



,(94) 



where Vd{k) is given by Eqs. ()76l77|) . For temperatures 

firmn = " mm[J A{k) + Qvd{k), JA{k) - 2Qw<j(fc)]. (95) 

The Fourier transformed dipolar interaction kernel is pos- 
itive definite, Vd{k) > 0. An unscreened dipolar interac- 
tion leads to a Vd{k) that diverges (tends to negative in- 
finity) at its minimum at fc = 0. In the presence of both 
upper and lower distance cutoffs (see, e.g., Eq.® for a 
lower cutoff) on the dipolar interaction, the minimum of 
Vd{k) attains a finite value and the system has a finite 
critical temperature. 



Examining Eq. ((94|) , we see that the introduction of the 
angular dependence in the dipolar interaction changes 
the results that would be obtained if the angular depen- 
dence were not included in a dramatic way. 

(i) New correlation and modulation lengths arise from 
the second term in Ea. ([M)) . 

(ii) At low temperatures, the second term in Eq. ((94|) 
becomes dominant as its poles have a smaller real part 
(and thus a larger correlation length) relative to the first 
term in Eq. fM]) that appears for an isotropic dipole- 
dipole interactions. 



D. Dzyaloshinsky- Moriya Interactions 

As another example of a system with interactions hav- 
ing non-trivial directional dependence, we consider a sys- 
tem of three component spins with the Dzyaloshinsky- 
Moriya interaction |351] present along with the ferromag- 
netic interaction and a long range interaction, 

H = -jJ2 S{x) ■Siy)+Y,D- [Six) x S{y)] 
{S,y) {S,y) 

+Qj2VLi\x-y\)S{x)-S{y). (96) 

We diagonalize this interaction kernel to obtain a Hamil- 
tonian of the form, 



^ = EE^a(2^)K(x,y)^a(y). 



(97) 



x,y a 



The Sa's are linear combinations of the components of 
S. In a large n approximation, the two point correlator 
is given by 



G(f) = kgT 



(Pk 



JA{k)+QvL{k)+ n 



2(JA(fc) + QwL(fc)+/i) 



(JA(fc) + QvL{k) + + {Dl + DI + Z?i)[A(fc)]2 



The presence of the Dzyaloshinsky-Moriya interaction 
does not alter the original poles and hence does not 
change the original lengthscales of the system. However, 
additional lengthscales arise due to the second term in 
Eq.i 



VIII. CONCLUSIONS 

In conclusion, we 

(1) showed the existence of certain crossover tempera- 
tures at which the lengthscales of the system change char- 
acter. As the long range interaction strength is turned 
off, the crossover temperature tends to the critical tem- 
perature of the unfrustrated ferromagnetic system. The 
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divergence of the modulation length at the crossover tem- 
perature defines a new critical exponent vi^. We proved 
that, quite universally, vl = 1/2. 

(2) studied the evolution of the ground state modula- 
tion lengths in frustrated Ising systems. 

(3) discovered a divergence of correlation lengths in 
the high temperature limit of long range large n sys- 
tems. This divergence is replaced by a saturation when 
the long range interactions have are screened. Notwith- 
standing this divergence, taking note of the prefactors 
associated with these correlation length, the amplitude 
of the correlations themselves decreases with increasing 
temperatures. 

(4) investigated, in large n theories, the evolution of 
modulation and correlation lengths as a function of tem- 
perature in different classes of systems. 

(5) proved that, in large n theories, the combined sum 
of the number of correlation and the number of modu- 
lation lengths is conserved. We have also showed that 
there exists a diverging modulation length at high tem- 
peratures for systems with long range interactions. 

(6) studied three dimensional dipolar systems. We 
found that the full dipolar interactions with angular de- 
pendence included, changes the ground state of the sys- 
tem and also adds new lengthscales. 
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with Lincoln Chayes, Daniel Kivelson, Steven Kivelson, 
Joseph Rudnick, Gilles Tarjus, and Peter Wolynes. We 
also thank the Center for Materials Innovation - Wash- 
ington University in St. Louis for providing partial fi- 
nancial support for this project. 



APPENDIX A: TRANSFER MATRIX IN THE 
ONE-DIMENSIONAL SYSTEM WITH ISING 
SPINS 



positive and is non-degenerate. The secular equation, 
det(T — A/) = is a polynomial in A with real coeffi- 
cients. Thus, two possibilities need to be examined: real 
roots, and, complex conjugate pairs of roots. Real eigen- 
values \p give terms with correlation length. 



(A2) 



Complex conjugate eigenvalues, \q and A* correspond to 
the same correlation length and modulation length, given 

by, 



^ = In 



\\\ 



L, 



27r 



tan 



ReXn 



(A3) 
(A4) 



Thus, the total number of correlation and modulation 
lengths is the order of the polynomial in A in the sec- 
ular equation, or simply the dimension of the transfer 
matrix - 0(2^). Similar to our conclusions for the high 
dimensional continuous spin systems, this number does 
not vary with temperature. For q state Potts type spins, 
replicating the above arguments mutatis mutandis, we 
find that the total number of correlation and modulation 
lengths is 0{q^). Similarly, for such a system placed on 
a d dimensional slab of finite extent in, at least, {d — 1) 
directions along which it has a length of order 0{l) > R, 
there will be 0{q^) transfer matrix eigenvalues and thus 
an identical number for the sum of the number of mod- 
ulation lengths with the number of correlation lengths. 

The eigenvalues change from being complex below cer- 
tain crossover temperatures to being purely real above. 
These temperatures form the "disorder line" . 



Thus far, we focused primarily on high dimensional 
continuous spin systems. For completeness, we review 
and illustrate how some similar conclusions can be drawn 
for one dimensional Ising systems with finite ranged in- 
teractions and briefly discuss trivial generalizations. In 
particular, we show how the sum of the number of modu- 
lation and number of correlation lengths does not change 
as the temperature is varied. In Section ljVI C 3|) . we il- 
lustrated how this arises for general large n systems. 

For interactions of range i? in a one dimensional Ising 
spin chain, the transfer matrix, T is of dimension AI = 
0{2^). The correlation function for large system size, 
takes the form 



k=l 



Ar 



(Al) 



where A.^s are the eigenvalues of the transfer matrix. 
Since the characteristic equation has real non-negative 
coefficients, from Perron-Frobenius theorem, Aq is real. 



appendix b: detailed expressions for 
sld for different orders p > 2 at which 

the long range interaction kernel 
has its first non- vanishing derivative 

If the lowest order non- vanishing derivative of VL{k) 
at the minimizing wavenumber fco (see Ea. (l37|) ). of order 
higher than two is that with p ~ 3 (the general case) 
then, in the large n limit, the change in the modulation 
length at temperatures T > Tc about its value aXT = Tc 
of Eq. (|i7)) is given by 



27r Qv'l\kQ)5^i 
^0 12(J 



5L 



D 



(Bl) 



We employ Eq. (|Bl[) in our analysis in Section ljVIip . If 
the lowest order derivatives are of order p = 4 or 5 then. 



SLd = 



2tt Qv^-^\ko){SfiY 
kl 24o(ji^!£!M^3 



(B2) 
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Similarly, for p = 6 or 7, 

^0 10080(J+ ^^^ijr^)^ 

and so on. 

APPENDIX C: tJ-{T) FOR THE SCREENED 
COULOMB FERROMAGNET 

In three dimensions, with A = 27r/a an ultra-violet 
cutoff with a the lattice unit length, at high temper- 
atures, T > T*, this leads to the following implicit 
equation for fi{T) in the case of the screened Coulomb 
ferromagnet of Eg. ([75)1 . 



1 _ A ^2 

X , = tan ( , 

-^!^i^£±^±^ tan-(^^)) (CI) 

Here, we employed the shorthand p = \/{p.— ~ 4Q- 
This parameter p vanishes at the crossover temperature 
T* at which a divergent modulation length makes an ap- 
pearance, p{T = T*) = 0. At low temperatures, T <T*, 
p becomes imaginary and an analytical crossover occurs 
to another real functional form. 
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